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It is theoretically known that a pair of phase conjugating surfaces can function as a perfect lens, 
focusing propagating waves and enhancing evanescent waves. However, the known experimental 
approaches based on thin sheets of nonlinear materials cannot fully realize the required phase con- 
jugation boundary condition. In this paper we show that the ideal phase conjugating surface is in 
principle physically realizable and investigate the necessary properties of nonlinear and nonrecipro- 
cal particles which can be used to build a perfect lens system. The physical principle of the lens 
operation is discussed in detail and directions of possible experimental realizations are outlined. 

PACS numbers: 42.65.Hw, 78.67.Pt, 81.05.Xj 

I. INTRODUCTION 

The perfect lens 1 is a device which focuses the field of a point source into a point, that is, the perfect lens focuses both 
propagating and evanescent fields. It is known 1 that a planar slab of the ideal Veselago medium 2 with the relative 
permittivity and permeability both equal to —1 has the perfect-lens properties because of the negative refraction 
phenomenon and the excitation of coupled surface plasmon-polaritons at the slab interfaces. Practical realization of 
such double-negative (DNG) materials is, however, a significant challenge, especially at optical frequencies, and, clearly, 
any realization will suffer from some imperfections. For example, metal-dielectric metamaterials have relatively high 
ohmic losses that are responsible for nonvanishing imaginary parts in the constitutive parameters of such volumetric 
artificial media. 

There is, however, a possibility for a different realization of a perfect lens that does not require any volumetric 
metamaterials. Indeed, the physical effects necessary for perfect lensing happen at the lens interfaces, and not within 
the metamaterial volume. Therefore, if one realizes a metasurface at which the incident waves refract negatively, 
then a parallel pair of such planar sheets will mimic the operation of the Veselago lens for the propagating plane 
waves. Moreover, if this metasurface supports surface modes (surface plasmon-polaritons) within a wide range of the 
tangential propagation factors k t > fco (where ko is the free space wavenumber), then also the impinging evanescent 
plane waves will interact resonantly with the sheets and will be tunneled through the lens with an enhanced amplitude, 
due to the electromagnetic coupling between the surface states excited on the sheets. Such subwavelength imaging 
with linear plasmon-polariton resonant grids was theoretically predicted and confirmed experimentally in a number 
of works. 3-7 

In 2003, we showed 8 that two parallel sheets with phase conjugating boundary conditions for tangential fields on 
the two sides of the sheets 

E t+ = E*_, H t+ = Hj_ (1) 

have the necessary properties of the perfect lens outlined above. In these conditions that are written for the complex 
amplitudes of the time-harmonic fields (symbol * denotes complex conjugation operation) the indices ± indicate the 
field values on the two sides of an infinitely thin phase conjugating sheet. 

Obviously, boundary conditions (1) cannot be realized using linear materials, and in the same paper 8 the use of 
three- wave mixing in a nonlinear layer was proposed as an approach to realization of this effect. Phase conjugation and 
"time-reversal" devices were studied also earlier for other applications. It is interesting that in the same year (2003) 
an experimental microwave realization of a phase conjugating layer was published 9 independently from our work. 8 
Later, the concept of perfect lensing based on two nonlinear sheets was studied theoretically in Ref. 10. Alternative 
experimental realizations of nonlinear negative refraction effect were published in Refs. 11,12. 

However, in known devices based on antenna arrays with mixers, 9 or sheets of nonlinear dielectrics, 10 or arrays of 
only electric or only magnetic particles with nonlinear insertions, 12 the phase conjugated ( "time- reversed" ) products 
create waves propagating symmetrically to the both sides of the sheet, i.e., there appears a retrodirected wave 
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FIG. 1: An ideal Veselago lens: A planar slab of a double-negative (DNG) material with the medium parameters e = — eo and 
/i = — /io in free space. 



propagating back to the source. While the perfect lens operation can be theoretically approached even in such 
systems if the amplitudes of the nonlinear products tend to infinity in the assumption of nonphysical infinitely strong 
external pumping, 10 the ideal phase conjugating boundary conditions (6) cannot be realized within this scenario. 

In this paper we discuss the physical meaning of the ideal complex-conjugation boundary conditions (1) and outline 
possible approaches for realization of such surfaces, which would potentially lead to creation of super-resolution lenses. 
The paper is organized as follows. In Section II we consider physical processes taking place at a phase conjugating 
boundary and demonstrate that such boundary may be equivalently represented with pairs of electric and magnetic 
surface currents reacting (nonlinear ly) to the applied magnetic and electric fields. In Section III a realization of the 
phase conjugating boundary with an array of bi-anisotropic inclusions is proposed and studied and the necessary 
conditions on nonlinear susceptibilities of the inclusions are established. In Section IV a possible microwave design of 
such inclusions is proposed. 



II. THE PHYSICAL MEANING OF THE COMPLEX-CONJUGATING BOUNDARY CONDITIONS 



A. Complex-conjugating boundary and perfect lens 

Let us start from outlining the idea from our paper in Ref. 8. Consider an ideal Veselago lens depicted in Figure 1. 
Let the relative permittivity and permeability of the medium surrounding the lens be equal to I and the relative 
parameters of the lens material to —1 at the working frequency to, respectively. At the lens interfaces the tangential 
components of the fields satisfy the usual Maxwellian continuity conditions. One may notice that in this system 
the only difference between the time-harmonic [the time dependence is of the form exp(+jut)] field equations in the 
Veselago slab (region 2) 

VxE = jw^oH, V x H = -jweoE (2) 

and the analogous equations in the free-space regions is the sign in front of the imaginary unit. A substitution 

E(oid) , H (old) => CE( ncw) , CH( ncw) (3) 

(C is an arbitrary constant; here and thereafter * denotes the complex conjugation) into the field equations in region 
2 reformulates the problem in terms of the new field vectors in which the field equations become the same in all three 
regions: 

V x E = -jufiotl, VxH = jwe E, (4) 

that are simply the Maxwell equations in free space. The boundary conditions on the two interfaces, however, are no 
more the standard continuity conditions, but they involve complex conjugation: 

E t (1 , 3) = C E t (2) : H t (1|3 ) = CH t * (2) . (5) 

The constant C describes the "transformation efficiency" of the nonlinear surface which transforms fields into the 
complex-conjugate state. In the known experimental realizations of phase conjugation in electrically thin layers (e.g., 
Ref. 12), the efficiency has been rather small. However, by choosing a small value of C in (5), we arrive at a structure 
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with asymmetric properties with respect to the two sides of the surface. Indeed, complex conjugating and dividing 
(5) by C* , we see that in this case the weak fields inside the lens should be enhanced by the surface in the same rate 
as they are suppressed when the surface is excited from outside. For this reason we will concentrate on the simplest 
choice of C = 1, as in Ref. 8, described by the boundary conditions 

E t(i, 3) = E t (2)5 H t(i, 3 ) = H t* (2) - (6) 

In this case the complex-conjugating surface has symmetric properties with respect to its two sides, and for the ideal 
lens operation the amplitudes of the field should not change across the sheets. 

Now it becomes evident that the problem involving an ideal Veselago slab is mathematically equivalent to the 
problem dealing with a pair of conjugating surfaces in free space, provided that the field sources are outside of region 
2. Therefore, in the latter system the field solutions are the same as in the Veselago slab, and because of this the 
physical phenomena taking place at the interfaces of region 2 are also the same: The propagating plane waves are 
refracted negatively at the interfaces, and the evanescent modes are enhanced due to the excitation of coupled surface 
plasmon-polariton pairs. In this regard, a pair of phase conjugating planes is indistinguishable from a perfect lens 
proposed by Pendry. 1 In what follows we concentrate on physical properties of such phase conjugating sheets. 



B. Plane-wave propagation across the phase conjugating sheet 



Let us consider a single phase conjugating surface located at z = 0. We decompose the tangential electric and 
magnetic fields into plane waves at both sides of the surface: 



E t (x,y) = J— [ [ Et (k x ,k y ) e -J(fc.*+M) dk x dk y , 

z=±o y^i J J z=±o 

H t (x,2/) = — ^ f fu t (k x ,k y ) e-^ k * x+k ^ dk x dk y . 

z=±0 (27Tj z J J z=±0 

It is easy to see that the boundary conditions (6) require that the plane-wave components satisfy 



Ht(fc x , ky) 



z=+0 



z=+0 



— E t ( k x , k y ) 
= H t (— k x , —k y ) 



z=-0 



(7) 

(8) 

(9) 
(10) 



From these relations we immediately realize that the propagating modes refract negatively at the conjugating interface 
due to the change in the sign of the tangential component of the wave vector k t = fc^xo + k y yo- This is very different 
from the case of the same refraction at an interface with a double-negative medium where the normal component of 
the wave vector changes sign. 

What about the evanescent modes? It can be shown that they are at resonance with the phase conjugating surface 
so that the strong (theoretically infinite) reflection takes place at the surface, although the field transformation in the 
sheets (the complex conjugate operation) does not amplify the fields. The details are given in Ref. 8; here we will just 
try to convince the reader with the following simple physical argument. 

Let us consider a single plane- wave component interacting with a single phase conjugating surface in free space. 
The plane wave is incident from the half space z < 0. On both sides of the phase conjugating sheet the tangential 
fields satisfy the usual relation between the fields in a free-space plane wave (the following equations are valid for TM 
or TE waves separately): 



=±o 



-Z T m,te zo x H t (k x ,k y ) 



E t (k x ,k y ) 

Here the free-space wave impedances for TM- and TE-polarized waves read 

1 



z=±0 



Ztm = ?jo y 1 — £2 ' 



Zte — Vo 



1 



k 2 



(11) 



(12) 



where ko = ui/c is the free-space wavenumbcr, rjo = yTXo/eo, an d &t = + by- Bvrt on the other hand the tangential 
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fields satisfy also the boundary condition (6). Together with the impedance relation (11) this leads to 



Et(k x ,ky) 



E* t (-k x ,-k y ) 

Z = + 

"■^TM.TE z o x H t ( — k x , —k y ) 



^TM.TE z x Ht(k X7 ky) 



z=+0 



(13) 



Thus, a wave incident on the sheet from one side "sees" the surface impedance which equals to the complex conjugate 
of the wave impedance in free space. The reflection coefficient reads 



R 



TM.TE 



^TM.TE — ^TM,TE 
^TM.TE + ^TM.TE 



(14) 



If the wave is a propagating wave, that is, |k t | = k t = ^Jk% + k^ < fc = uj/c, then its wave impedance is a 

real number, and the reflection coefficients equal zero. This proves that the propagating modes experience negative 
refraction without any reflection. 



If the wave is evanescent, that is, fc t = yk^ + ky > fco, the wave impedances (11) are purely imaginary and the 

reflection coefficient is infinite. The transmission coefficient is also infinite due to the boundary condition (5). This 
is obviously a resonant condition which can be also understood as a condition for existence of a surface mode. 8 Note 
that this resonant condition holds for all values of the tangential wave number k t > fco, which is the condition for 
perfect lensing of all evanescent field components. 

The physical reason of such a resonance is rather simple. Consider, for instance, the waves of TM polarization. 
The characteristic impedance of an evanescent TM wave (12) can be written as 2tm = — jVo®/ko, where a is the 
decay factor: a = \/k\ — k$. We conclude that a TM evanescent wave has capacitive characteristic impedance. Due 
to the conjugating interface, the same impedance of the matching TM wave behind the interface (at z = +0) is seen 
as inductive in front of the interface (at z = — 0), which reminds of a "connection" of two reactances of the opposite 
character in a parallel oscillatory circuit. Probably, A. Alii and N. Engheta 13 were the first to identify and explain 
the resonance of the same nature that happens at the border of a double-positive and a double-negative material. 

Therefore, when an incident evanescent wave excites a conjugating surface, it resonantly excites a surface mode 
that matches its transversal propagation factor and this results in very strong (theoretically infinite) reflected and 
transmitted waves at the surface. The strong reflection was also found to be the key to sub-wavelength imaging in a 
pair of parametrically pumped nonlinear nonmagnetic sheets studied in Ref. 10 The above discussion shows that in 
a metasurface realizing the boundary conditions (6) the reflection and transmission coefficients for evanescent modes 
tend to infinity due to a high-quality resonance (theoretically, with an infinite quality factor), while, in Ref. 10, the 
surface itself parametrically amplifies the fields. 



C. Equivalent surface currents on the phase conjugating sheet 



One may notice that the boundary conditions (6) imply discontinuity of both tangential electric and magnetic fields 
across the phase conjugating plane. The jumps of the fields can be expressed as follows: 



H t (x,2/) 



z=+0 



z=+0 



z=-0 



z=-0 



-2jIm[E t (z,t/)] 
-2jlm[U t (x,y)] 



z=-0 



z=-0 



These jumps are related to the equivalent magnetic and electric surface currents that exist on the surface: 



z x J m (x, y) = E t (x, y) 



z=+0 

-2jIm[E t (x,2/)] 



-z x 3 e (x,y) = H t (x,y) 



z=+0 



Et(x,y) 

z=-0 

Ht(£,y) 



z = -0 

2jIm[E t (a;,*/)] 

z=-0 



= -2jlm[n t (x,y)} =2jTm\H t (x,y)] 

z= — 



z=+0 



z=+0 



(15) 
(16) 



(17) 



(18) 
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Let us stress already at this point that the relations (17)-(18) are the only physical conditions one has to satisfy 
in a subwavelength imaging device based on phase conjugating sheets. Nothing more is required! Essentially, these 
relations tell us that it is enough in practice to realize a metasurface that reacts with certain magnetic and electric 
surface currents to the imaginary parts of the tangential electric and magnetic fields at a given side of the surface. 
This also shows that there is no need for extra-strong pumping or super-efficient nonlinear conversion as in Rcf. 10. 
Moreover, from (17)-(18) we observe that the induced electric current should be proportional to the magnetic field on 
the sheet (more precisely, to its imaginary part). Likewise, the induced magnetic current is proportional to the electric 
field. This is quite different from the known approaches based on layers of nonlinear dielectrics or magnetics. 10 ' 12 

Let us decompose each of these surface currents into a sum of two currents: J c = Jq 1 ' + ji 2 \ J m = Jm + Jm\ 
where 



-z x J« = -H t , z x J« = -E t 

z= — 



(19) 

z= — 



-zo 



x J< 2 > = H t 



, z x j£> = E t 

z=+a 



(20) 

z=+0 



One can see that the pair of the surface currents Je , Jm is essentially an equivalent Huygens source defined at the 
plane z = — 0. In the half-space z > this source produces the field which is the negate of the field that the external 
sources located at z < induce in the half-space z > (the negation is due to the minus signs in the right-hand 

side of (19)). Thus, the physical role of the currents Jm' when concerned with the half-space z > is to cancel 

(2) (2) 

the field incident from the other half-space. The same holds for the other pair of currents Jo , J m when concerned 
with the half-space z < 0. These currents form a Huygens source defined at z — +0 plane. In the half-space z < 
they cancel (pay attention to the direction of the normal!) the field produced by the external sources from the z > 
half-space. 

From the other hand, the pair of currents ji 1 ^ , Jm plays another role when concerned with the half-space z < 0. 
Indeed, using the boundary conditions (6) we write 



-zo x 



J« = H 



j!! } - -e: 

z=+a 



z x „ m - 



(21) 

z=+0 



from which it is evident that these currents can be identified also as an equivalent source located at the plane z = +0 

(2) (2) 

that produces at z < the conjugated field of the sources located in the half-space z > 0. Respectively, J e , J m 
produce the conjugated field in the region z > 0. 

To summarize, we have identified the following roles of the currents: 

• The pair Je 1 ', Jm"* cancels the field incident from z < to the half-space z > and creates the conjugated field 
of the sources from the half-space z > in the half-space z < 0; 

(2) (2) 

• The pair J , JV cancels the field incident from z > to the half-space z < and creates the conjugated field 
of the sources from the half-space z < in the half-space z > 0. 



III. PHASE CONJUGATING SURFACE AS AN ARRAY OF BI-ANISOTROPIC NONLINEAR 

INCLUSIONS 

From the above results we see that the ideal phase conjugating surface should respond to the fields with both electric 
and magnetic polarization. Dependence of the induced electric current on the magnetic field and vice versa suggests 
that the structure should have some magneto-electric coupling. In this section we investigate if it is possible to realize 
the ideal phase conjugating boundary conditions (6) with a planar array of nonlinear bi-anisotropic particles. 



A. General requirements on susceptibilities of inclusions 

Let us first find out how the total induced electric and magnetic surface current densities depend on the incident 
electric and magnetic fields in the array plane. To do that, we consider an isolated phase conjugating surface in the 
field of a single TM (or TE) polarized plane electromagnetic wave (propagating or evanescent). Taking into account 
the conjugating boundary condition (6) we may formally write the total tangential electric and magnetic fields on 
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both sides of the surface as 





= (l + R T M.Tv)K nc (x,y) 

z = — Q 


2=0 




Et(x,y) 


= (l + RTM.TE)(n nC (x,y 

z=+0 


))* 




H t (a;,y) 


= (l-i?TM.T E )Hi nc (x,y) 

z = — 






H t (a!,y) 


= (l-VTE)(Hi nC (^))' 

z=+0 


z=0 



where the reflection coefficients -Rtm.te are given by (14), from which we notice that -R TM TE = - 
expressions hold for both propagating and evanescent waves incident from the half space z < 0. 
Therefore, from (17)-(18) and (22)~(25) wc obtain 



(22) 
(23) 
(24) 
(25) 

-Rtm te- The above 



-zo x J m (x,y) 
z x J e {x,y) 



2jIm(Ej 
2jIm(H 



(z,2/)) 

\x,y)) 



z=0 



z=0 



2i?TM,TE R°(E" 

- 2_Rtm,te Re(H 



(x,y)) 

z {x,y)) 



z=0 



z=0 



(26) 
(27) 



The addends on the right-hand side of (26)-(27) that are proportional to the imaginary part of the incident field 
are relevant for the propagating waves (for these waves -Rtm.te = 0), while for the evanescent waves the addends 
proportional to the real part of the field are of the most importance, because -Rtm.te | — > °o for these modes. It is 
instructive to compare these observations with the discussion in Sec. II C. From (17)— (18) it follows that in terms 
of the total tangential fields at a given side of the phase conjugating sheet, the conditions for both propagating and 
evanescent waves are the same: Eqs. (17)-(18) do not distinguish these waves. Physically, this is because the locations 
where one must measure the fields and where one must create the surface currents are not at the same point, if one 
wants to approach a design directly suggested by these equations. Indeed, the mathematical form of Eqs. (17)— (18) 
demands that the field values must be taken at a point slightly displaced off the surface z = 0. On the contrary, in 
an array of particles (considered here as point objects) the equivalent surface currents depend only on the fields in 
the array plane and, as we will see soon, the required reaction to this field happens to be different for the two types 
of waves. 

A related observation is that in a realistic structure, e.g., an array of nonlinear polarizable bi-anisotropic inclusions, 
it is not the incident field, but the local field E[ oc , H[ oc , that excites each and every inclusion in the structure. The 
latter has a contribution from the secondary field of the induced currents. We may therefore write for the signals at 
the frequency of the incident wave: 



(x,y) 



Hi° c (x,y) 



-pane 
jjinc 



(x,y) + P ee - J e (x,y), 
{x,y)+]3 mm • J m (x,y), 



(28) 
(29) 



where f3 cc mm are the so-called interaction dyadics. For arbitrary distributed currents these dyadics are understood 
as operators acting on the currents. However, for the following it is enough to consider the currents of the form 
Je,m(a;, y) — Jk™ ex P(~ jk t • r) + J^™ exp(jk t • r). In this case, in order for (28)-(29) to hold in a simple dyadic sense, 



the interaction dyadics must satisfy /3 ee mm = /3 e e,mm(kt) = /?cc.mm(~k t ), i.e., the lattice (not the particles!) must 
have a center of symmetry. There are no cross-terms in (28)-(29) because the tangential magnetic (electric) field of 
an array of tangential electric (magnetic) dipoles vanishes in the plane of the array. 

Additionally, the induced electric and magnetic currents must be sensitive to the phase of the external field, because 
the conjugating boundary reacts differently to the real and imaginary parts of the tangential electric and magnetic 
fields. Therefore, the inclusions must react differently to the corresponding components of the local fields. Based on 
the above discussion we may write 



J c = a° c ■ Re(E. 



loc 

J m = Ce • Re(Ei° c ) + 02 ■ Im(El° c ) + 1 



) + Me ' Im ( E 



loc^ 



Re(U l n+0Z 



:im 
mm 



Im(Hl oc ), 



(30) 
(31) 



where Q!gg'™ m mc cm are the dyadic polarizabilities to the real and imaginary components of the local fields. 
Substituting (31)-(30) into (28)-(29) we obtain 



(I t -(3 cc -of c c c )-MK c )+3(It~P c 



— im > 



Im(Ei oc ) 



Rc(Hi oc ) - jfj m ■ <C • Im(Hl° c ), (32) 
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H™= = (h - /3 mm • <Cn) • Rc(Hi oc ) + j(I t - (3 mm ■ a£n) • Im(Ht oc ) 

- L m • Co • Re(E[° c ) - mm ■ 1™ • Im(Et oc ), (33) 

where I t is the unit dyadic in the plane of the array. Respectively, 

Re(Er) = Re(7 t - ^ cc • % e e ) ■ Rc(E;° c ) + Im(^ ee • !Q • Im(Ei° c ) 

- Re(^ ee • Cm) • Re(Hi° c ) + Im(^ cc • 5^) ■ Im(Hl° c ), (34) 

Im(Ei nc ) = - Im(^ ee • C) • Re(Ej oc ) + Re(7 t - % c ■ !Q ■ MEi oc ) 

- Im(| ee • EZ) ■ Rc(H l n Re(% e ■ 5^) • Im(Hi° c ), (35) 

Re(Hl nc ) = Re(7 t - ^ mm • Cm) ■ Rc(Hi° c ) + Im(^ mm • i^J • Im(Hi° c ) 

- Re(^ mm • l^c) • Re(El oc ) + Im(^ mm • 5^) • Im(Et° c ), (36) 



Im(Hi nc ) 



- Im(/3 mm • <£j • Re(H[ oc ) + Re(J t - /3 mm • <C m ) • Im(Ht° c ) 

-Im(L m -Ce)-Re(Et oc )-Re(^ 



"3 =im 



l„;;).Im(Et oc ). (37) 



These expressions can be substituted into (26)~(27) from which one obtains a set of dyadic relations for the polar- 
izabilities assuming that the four components of the local fields Re(E\ oc (x , y j) , Re(H[ oc (a;, yj), Im(E[ oc (a;, yj), and 
lm(H l ° c (x,y)) are independent. Doing so we obtain the following relations: 
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=im 
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(38) 
(39) 
(40) 
(41) 
(42) 
(43) 
(44) 
(45) 



where R is the dyadic reflection coefficient defined in terms of -Rtm,te as 

= znZn^ktkt ktkt 

it — JlTE P5 ' -KTM " 



(46) 



(for the definition of the double cross product and other dyadic algebra rules see, e.g., Ref. 14). One may notice that 
because R in the above relations is either zero or purely imaginary: (R)* = —R, it follows that Re(a™ mmemme ) = 
Im(a™ mm em mc ) = 0. Therefore, the equations (38)-(45) can be also written as 
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Let us consider first the propagating part of the spectrum. For such waves, R = 0, and the above relations simplify. 

Also, we can write /3 ee = 770/8, and /3 mm = tJq 1 /3, where (3 is the dimensionless interaction dyadic which, by duality, 
is the same for the electric and magnetic currents as they are due to the electric and magnetic dipole moments that 
belong to the same particles in the array. The solution of the system of dyadic equations (47)-(54) in the case of 

R = is 
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cc,mm,cm,mc 



= o, 
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ma£ = Vo ^nTm = 4z o x Re(/3) 



(z x Re(/3)) 
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7 t + 4(z x Re®)' 
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It can be shown 15 ' 16 that the real part of the interaction dyadic (3 for a planar array verifies 

1 z z *k t k t . 



Re(/3) 



2cos6> 



m a., cos 6* k t k t fcgAoy 
"1 fcT + _ 6^ *' 



(55) 
(56) 

(57) 



(58) 



where 9 is the angle of incidence: cos 9 = y/l — k^/k^, and A is the unit cell area. This result holds for arrays with 
arbitrary unit cell geometries, provided that the arrays do not produce higher-order diffraction lobes. 

It is quite interesting that the imaginary part of the interaction constant that contains the information about 
the microstructure of the array has completely disappeared from the above solution. The imaginary part of the 
interaction constant does not contribute in this case because the induced currents J and J m are always imaginary, 

and the respective additions to the interaction field j Im(/3 ee ) • J e and j'Im(/3 mm ) • J m are real- valued, to which the 
particles do not react. Thus, the interaction of the particles in the array is irrelevant in the considered case, and each 
particle radiates effectively as in free space. 

We may substitute (58) into (56)-(57), taking into account that 



(z xRc03)) 2 = 



1 

'I 



klA 
6ir 



+ 



kjA Q 
12tt 



cos 9 + 



cos 9 



and obtain 



— im 



-a 



6tt / „ 1 

cos 9 + 



cos 9 



kjA Q 
3tt 



(z x I t ) 



-l=im 
^0 a mm 



127T 

k%A 



1 



cosf 



kjA 

cos 9 3tt 



(z z * Re(/S)). 



(59) 

(60) 
(61) 



For practical purposes, considering the phase conjugation of paraxial beams in dense arrays, we may approximate the 
above relations as 



3tt 

Wa 



— im _i=im 

Vaa cc =Vo "mm 



(z x I t ), 



3 '1. 



kM 



(62) 
(63) 



Because of the form of the relations (22)-(25), the obtained exact solutions (56)-(57) and their approximations (62)- 
(63) arc valid for arbitrary plane waves incident from the half space z < 0. The solution for the case of incidence from 

the half space z > is obtained by replacing z with — z in (56)-(57) and (62)-(63), which changes signs of a™ me . 
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From the above results we see that the particle must be "invisible" for the real- valued electric and magnetic fields, 
while the polarizabilities of the particle to the imaginary-valued fields must be such that the electric and magnetic 
currents form Huygens pairs that absorb the incident wave and produce the phase conjugated wave. 

= -i 

For the evanescent waves -Rtm,te| — > oo, therefore, it is convenient to multiply the equations (47)-(54) by R 
from the left. Then, in the limit R — > the following solution of the system (47)-(54) can be immediately found: 



— r r. _l=rc 



Im(/?) 



— im 



—ream 



0. 



(64) 
(65) 



The same solution can be also obtained with a more accurate treatment. Let us introduce the notations C = 
Rc(/3 00 ) +jR ■ Im(/3 ee ) and C m = Re(/3 mm ) — jR ■ Im(/3 mm ). Then, in these notations, we may, for example, write the 
solution for a no as 



a ee = 4 



It + 4(z x C m ) • (z x C e ) • (z x C m ) • (z x R) 



I t + -(z x C ) 1 • (z x C m ) 1 



(z x C ) 1 • (z x R) 



= C c 1 -i? + 0(i? 2 ). (66) 



Next, 



c/ • R 



ImCSee) 



jR 1 • Re(/3 co ) • (lm(/? ee )) 



(67) 



which leads to (64). 

From (64)-(65) we conclude that to conjugate the evanescent part of the spectrum the inclusions must be "invisible" 
to the imaginary part of the electric and magnetic fields. The inclusions do not have to be bi-anisotropic in this case. 
The particles are purely reactive and their reactance should compensate the reactance due to particle interactions, 
creating a resonant structure. 

From a physical point of view, condition (64) can be understood as a condition for a surface polariton resonance 
at the array surface. Indeed, for particles with the polarizabilities (64)-(65) we may write J c = a cc ■ Re(E loc ). From 
the other hand, E loc = E inc + /! cc • J c . Therefore, 



J = C • Rc(E inc + p cc ■ J c ) = cC • Rc(E inc ) + jaZ ■ Im(/? ee ) • J c 
because both a™ and J e are purely imaginary. From here 



(68) 



/t-jC-M/U •< C c -Re(E inc ), 



(69) 



and we see that J c — > oo when condition (64) is fulfilled. A similar resonance is responsible for the enhancement of 
the evanescent waves in a pair of linear plasmon-polariton resonant grids studied in previous works. 3-6 



B. Electromagnetic properties of the particles forming phase conjugating sheets 

Although the principle of operation of the field-conjugating perfect lens is the nonlinear operation of complex 
conjugation of electromagnetic fields, it is interesting to observe that the particles which perform this operation arc 
characterized by linear polarizabilities with respect to the real or imaginary parts of the fields. The nonlinear nature 
of the particles is thus only in their selective sensitivity to either real or imaginary parts of the complex amplitude of 
the local fields. 

Considering the particle response to the real or imaginary field components separately, we may apply the theory 
of usual linear bi-anisotropic particles. For the particles which react to the imaginary parts of the field (the particles 
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excited by the propagating part of the spatial spectrum), we rewrite relations (30), (31), and (62) in terms of the 
induced electric and magnetic dipole moments of individual particles p c ,m an d the local fields: 



-j—^ j j Im(Ei oc ) + Vo {j—^ j zo x j Im(H[° c ) 



- a cc ■ j Im(E l ° c ) + Voa cm ■ j Im(Ht° c ), (70) 



= I mm • j Im(Hi oc ) + -I mc • j Im(Ei oc ) (71) 
Vo 

(the surface current densities are related to the dipole moments of individual particles as J e , m = jwp C;m /io). 

These relations show that the particles reacting to the propagating part of the spectrum are bi-anisotropic and 
nonreciprocal. The magnetoelectric coupling is due to nonrcciprocity only (no magnetoelectric coupling due to 
reciprocal spatial dispersion effects), because the coupling dyadics satisfy 

a cm = a mc . (72) 

Furthermore, because these dyadics are antisymmetric (a cm = — a em , a mc = — 5 me ), materials formed by particles of 
this type belong to the class of moving media. 17,18 

The polarizabilities of lossless bi-anisotropic particles satisfy the following conditions (e.g., Ref. 18): 

^cc — ^eei ^mm — ^mm) ^cm ^mci (^^) 

where f denotes the Hermitian conjugation operation. Obviously, the inclusions with the polarizabilities (70) and (71) 
have the opposite property of being purely passive or active (there is no stored electromagnetic energy in their near 
fields), because they satisfy the opposite conditions: 

O-cc a cc> ^mm — ^mm? fl em — %e' (^) 

The power extracted from the local fields by one pair of the particles reads 

P = ^Re{J* • E loc + • H loc }A . (75) 

Assuming paraxial wave propagation, we can use the polarizability expressions (62) and (63) and assume that the 
electric and magnetic local fields are related by the free-space wave impedance 770- In this approximation we find, for 
the case of plane wave incidence from the half space z < 0, 

p = \%^ {EXoc)]2 = Wo llm(Eloc)]2 = 6 ~W [im(i/loc)]2 - (76) 

As is clear from this result, each of the polarizability components in (62) and (63) brings equal contributions to 
the extracted power. Noting that the induced dipole moments of ideal absorption-free dipole scatters read (at the 
resonance) 



t^Ioc „ /' .6717x0 \ „], « 



Poo = {-j-jf ) ' P-o = {-j-jjf ) H ( 77 ) 

we see that a pair of such ideal dipoles would extract from the fields exactly the same amount of power as our phase 
conjugating particles (when the complex amplitude of the local field is purely imaginary at the point of the particle). 
Thus, we can conclude that the particles described by (62) and (63) actually do not absorb power. They act as 
ideal absorption-free particles, which receive power from the incident field and re-radiate the same amount of power, 
creating phase conjugated waves of the same intensity as the incident propagating waves. 

It is easy to check that the same particles with the polarizabilities (70)-(71) do not react to the plane waves incident 
from the half space z > 0. Indeed, the sign of the magnetic field in these waves is opposite, therefore, the contributions 
due to E loc and H loc to the electric dipole moment p e and the magnetic dipole moment p m compensate each other in 
Eqs. (70)-(71), so that p e = p m = under such excitation. As was mentioned in Section III A, to conjugate the waves 
incident from the half space z > 0, one must change the signs of the magnetoelectric coupling terms in (70)-(71). 
Physically, this requires another array of inclusions with slightly different topology (more details in the next section) . 
Fortunately, the particles in the two arrays do not interact, so that in practice it is possible to combine the two types 
of particles in a single plane, for example, in a chess board-like structure. 
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IV. DESIGN OF THE PHASE CONJUGATING BI-ANISOTROPIC INCLUSIONS AT MICROWAVES 

A. Case of propagating waves 

To approach the design of nonlinear phase conjugating inclusions at microwaves one may start from the ideas behind 
the well-known omega particle. 19,20 An omega particle is a combination of a short dipole antenna and a small loop 
antenna. The particle is planar, so that both the dipole and the loop may be printed on a single side of a printed- 
circuit board. In the most common variant of this linear and reciprocal inclusion the dipole is directly connected 
to the loop. The nonlinearity and nonreciprocity, thus, can be achieved if one inserts a nonlinear and nonreciprocal 
four-pole network between the two antennas. 

To identify what kind of network one may need, let us hrst briefly analyze how the linear omega particle reacts 
to the local electric and magnetic fields. We consider the case when the particle lies in the xz-plane with the dipole 
antenna oriented along the x-axis. In this geometry, the dipole reacts to the x-component of the electric field, E l ° c , 
and the loop reacts to the y-component of the magnetic field, H l ° c . 

The electromotive force (EMF) induced by the local field in the dipole can be written as 

fdip = h dip E l ° c , (78) 

where h^ v is the effective height of the dipole antenna. For a short dipole of the total length 21, the effective height 
is hdip — I. Respectively, the EMF induced in the loop by the magnetic field reads 

£ioo P = -jojfioSH 1 ™, (79) 

where S — irr 2 is the area of the loop. Under a normal plane wave incidence, H x ° c = ±E 1 x oc /r t0 (the two signs are for 
the two possible directions of incidence), therefore, 

£ioo P - T.j(kvS)E l ° c = Tjh loop E l : c , (80) 

where we have introduced the effective height of a small loop antenna h\ oop — k S. One may notice that the EMF 
induced in the loop is in quadrature with respect to the EMF in the dipole. Therefore, when the two antennas are 
directly connected, these EMFs add up, but never fully compensate (or fully complement) each other. In fact, in the 
most interesting case when h\ oop = hdip = I, the total EMF induced in the antennas is 

£tot = (lTj)lE 1 : c . (81) 
Respectively, the induced current at the point where the dipole connects to the loop is 

di P ~~ ~7 T~7 ' ^ 2 > 

■^dip T -^loop 

where Zdip is the input impedance of the dipole and Z\ oop is the input impedance of the loop. The induced electric 
and magnetic dipole moments are proportional to this current: 

- ^^°\, (83) 



jiJ Oj(Zdip + ^loop) 

/ 1 -r- A\ j2 PIOC 

Pm,y = M0/dipS = V0 ) 7 +J ' X y (84) 
^(Zdip + Aoop) 

where we use the fact that S = l/ko if h\ oop = /idi P = I- When the particle is at resonance, Zdip + ^ioop = 2i? ra d + i?i OS s, 
nt)k 2 l 2 

where i? ra d = — is the radiation resistance of a short dipole antenna (when h\ oop = h^ v both antennas have 
the same radiation resistance) and i?i OS s corresponds to the ohmic loss in metal, which we may neglect. From these 
relations we see that the induced dipole moments in the linear omega particle are in quadrature. 

However, from (70)- (71) it follows, first, that in the nonlinear particle which we want to design the contributions 
due to E loc and H loc in the expressions for both dipole moments must be in phase (for the wave incident from z < 0), 
and, second, the induced electric and magnetic moments themselves must be also in phase. In a symbolic language, 
we may say that the four-pole network that we insert between the two antennas must act in such a way that Tj m 
(83)-(84) is replaced with ±1, and r/o in (84) with —jrjo- As will be seen in a few moments, this can be achieved with 
microwave nonlinear circuits known as balanced modulators (BM). 
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FIG. 2: (a) Equivalent circuit of an idealized balanced modulator (BM). (b) A network composed of two balanced modulators 
connected through a low-pass filter (LPF). 



A balanced modulator is a three-port device such that there are two ports that may both serve as input and output 
(the ports are exchangeable due to the symmetry of the circuit; that is one of the reasons why the circuit is called 
balanced modulator), and the third port to which a voltage from a local oscillator is applied. The function that the 
BM performs is a multiplication of the signal applied to one of its input ports and the signal of the local oscillator. 
One may represent an idealized BM with an equivalent circuit shown in Fig. 2(a). The controllable current sources 
in the circuit depend on the instantaneous voltages at the ports as follows: 

h{t) = K 12 u 2 {t)n {t), (85) 
i 2 (t) = K 2lUl {t)ui {t), (86) 

where ui (t) is the voltage at the local oscillator port and u\ j2 {t) and ii, 2 (t) are the voltages at the other two ports. 
We need power conserving BMs which do not absorb or store power that is delivered to ports 1 and 2, hence, 
U\{t)ii{t) + u 2 (t)i 2 (t) = 0. From here, K\ 2 = —K 2i . The input and output resistances of the BM shown in the figure 
with dashed lines are assumed to be very large and are not taken into account. 

Consider now the network depicted in Fig. 2(b). In this network we have connected two BMs through a low-pass 
filter (LPF). One may imagine this LPF as a n-type CLC-filter with an inductor in a series branch and two capacitors 
in the parallel branches. For us, however, the only thing that is important here is that this LPF freely passes through 
the direct-current (DC) component and blocks all high-frequency components (the DC path through the filter is shown 
with dashed lines). 

The whole network is designed to operate with signals at the frequency to = 2-kJq — 2n/T , therefore, we may 
represent the voltages Ui t2 (t) as 

ui, 2 (t) = U{ c 2 cos tot - U[ m 2 smut. (87) 

We apply the local oscillator signal at the frequency lo and the phase ip — tt/2 to the first BM: ui ,i = — sinLot, and 
another signal at the same frequency and ip = to the second BM: ui Q , 2 — cosujt. Therefore, we may write for the 
DC currents i[ 2= (here (. . .}t denotes averaging over a period): 

»i= = K 21 ((£/J°coswt- U[ m sin ujt)(- sin ujt)) To = ^K 21 U[ m , (88) 

*2= = K i2 ((U r 2 c cos Lot - U 2 m sintot) cosLot) To = ^K 12 U r 2 c . (89) 

But as is dictated by the topology of the network, i' 1= = —i' 2= , therefore, 

U[ m = t/ 2 rc . (90) 

Next, we express the high-frequency currents at the ports 1 and 2: 

h(t) = (K 12 u' l= )(- sinut) = -I^sinuot, (91) 
h(t) = (K 2 iu' 2= ) cos cot — P 2 C cos Lot, (92) 

where 7™ = K\ 2 vl 1= and I 2 e = K 2 \u' 1= . But again, from the topology of the network, the DC voltages satisfy 
u' 1= — u' 2= , therefore, 

ir = -i™. (93) 
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FIG. 3: The topology of the two types of phase conjugating bi-anisotropic inclusions for the phase conjugating surface operating 
at the frequency fo- Each particle consists of a short dipole antenna with the impedance Zdi P and a small loop antenna with 
the impedance Z\ oop . The two antennas are interconnected through a nonlinear and nonreciprocal network composed of two 
balanced modulators (BM) and a low-pass filter (LPF). The mixers are fed by a local oscillator with two signal outputs in 
quadrature ip — 0° and (p — 90° . See the main text for further explanations. 



Notice also that always I[ c = = 0. 

Thus, from (90) and (93) it is evident that the considered network operates, essentially, as a "connector" between 
the imaginary current and voltage at the first port, and the real current and voltage at the second port. The network 
also enforces zero real current in the first port and zero imaginary current in the second port. This is exactly what 
we need in the design of the phase conjugating particles, and the corresponding topologies including the antennas are 
shown in Fig. 3. 

In these designs we connect the dipole and loop antennas to the BM-based network discussed above. The electrical 
size of the circuit is negligible, and both antennas are excited by the same local field. Let us analyze the operation of 
the variant shown at the top of Fig. 3. First of all, it is evident that when placed in an arbitrary local field, the EMFs 
in both antennas are still given by (78) and (79). However, contrary to what happens in a linear omega particle, the 
real part of £dip and the imaginary part of £i 00 p will not be able to excite any current in the antennas, because these 
components are blocked by the BMs. Therefore, the relevant parts of the EMFs in the two antennas are (as above, 
we let h dip = /i loop = I) 

£™=llm(E l x ° c ), (94) 
££ op = m lIm(H l ; c ). (95) 

Analogously, the voltage drops on the reactive parts X,a p = Im(Zdi p ) and X\ oop = Im(Zi oop ) of the input impedances 
of the dipole and the loop have no effect as well, as they are in quadrature with respect to the current flowing through 
them. Thus, the only relevant part of the input impedance is the radiation resistance. 

From the topology of the network, 7 d ^ = —7™ and If® = I™, therefore, I™ p — I\£ op - Next, because of (90), the 
EMFs (94) and (95) and the rest of the equivalent circuits of the two antennas (only two i? ra d remain) appear to be 
essentially connected in series, hence, 



i[Im(^ c )+^I m(gj,° C )] 

2-Rrad 



= iio OP = L y v n ■ 06) 



The complex amplitudes of the electric and magnetic dipole moments, therefore, read 



p c , x = — = ——{Im(Eln + m Im(H^)}, (97) 

JUJ ^WKrad 

Pm,„ = MoCopS = [Im(i# c ) + 770 Im(i7'° c )], (98) 
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FIG. 4: A nonlinear single-port network used as a load for a short electric dipole. The network is composed of the same 
elements as in Fig. 2, with an additional operational amplifier working in the current-to-voltage conversion mode. 

which is a particular case of (70)-(71) for the considered polarization. 

As can be readily verified, the second design variant shown at the bottom of Fig. 3 has the opposite signs of the 
magnetoelectric interaction terms, and, thus, must be used to conjugate the plane waves incident from the half space 
z > 0. 



Case of evanescent waves 



To phase-conjugate the evanescent waves, the particles must react to the real parts of electric and magnetic fields, 
as follows from (64)-(65). In this case the inclusions are simple electric and magnetic dipoles without magneto-electric 
interaction. Therefore, as the base for our design we may choose the loaded dipole and loop antennas. In what follows, 
we consider in detail the case of a linear electric dipole oriented along the x-axis (the magnetic dipole along the y-axis 
may be considered in a dual manner). 

When a short loaded dipole is placed in an electric field, an EMF is induced in the dipole, with the value given 
by (78). Respectively, the current induced in the dipole is 



iE l : 



Zdip + -^load 



(99) 



where Z\ oa( i is the impedance of a bulk load connected to the dipole. The induced electric dipole moment of the 
loaded dipole reads 



Pe.x — 



loc 



= -3 



w(Zdip + Zload) 



It is easy to verify that if we chose the load so that 



^load = -Zdip + ] -^— Im(/3), 

^0 



(100) 



(101) 



then the condition for the surface electric current polarizability (64) in an array of such particles will be satisfied, 
with an exception that the linear particles will react to both real and imaginary parts of the electric field. 

With the use of BMs we may get rid of the reaction to the imaginary part of the field and also find a simple way 
to realize the necessary loading (101). Consider the network shown in Fig. 4. This network is similar to the ones we 
considered in Sec. IV A. It is composed of a pair of BMs and a LPF, with an additional element in the middle that is 
an operational amplifier working in the current-to-voltage conversion mode. We pump the first BM at the frequency 
/o with the phase (p = 0, and the second one with the phase if = ±7r/2 (it will be seen soon why we consider two 
signs of the phase). The output of the network is directly connected to its input, so that the circuit is essentially a 
single port device that may be used as an active load. To analyze the operation of the circuit we first notice that 
because the voltage at the input of the operational amplifier is negligibly small, u' 1= — 0, and, therefore, i\ — 0. Next, 
assuming that the voltage at the input of the circuit is u\ = U\ c coscot — U[ m sincjt, lo = 2i:fa = 2n/To, we obtain 



1 



A= = K 2i ((Ur coscot - [7™ sincjt) cosut) To = -K 21 U{ C 
Hence, the DC voltage at the output of the operational amplifier is 

«' 2 = - h=R = \k 21 rut, 



(102) 



(103) 
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FIG. 5: The topology of the two particles loaded with nonlinear circuits designed to operate with the evanescent part of the 
spatial spectrum. Top: a short electric dipole loaded with a nonlinear load. Bottom: a small magnetic loop loaded with a 
nonlinear load. 

where R is the resistor in the feedback loop of the operational amplifier. Therefore, 

i 2 (t) = tK 2 iu' 2= smut = T^K^RUl c sin cut = -I^sincut, (104) 

where J™ = iK^RUl /2. This current is the input current of the whole network, because the input current of the 
first BM equals zero. Thus, we have designed a circuit in which a real input voltage induces an imaginary current, 
i.e., the circuit behaves almost as a usual reactance (the plus sign in the expression for the current corresponds to the 
capacitance and the minus sign — to the inductance), with the difference that the loading circuit is not sensitive at 
all to the imaginary input voltage. This is exactly what we need to realize the nonlinear particles reacting only to the 
real part of the electric field. Indeed, to realize the required loading one has to choose the parameters of the circuit 
so that 

±\k\ x R = -X dip + ^ Im(/3). (105) 

It is interesting to note that there is no need to compensate the real part of the dipole impedance Z^ip (the radiation 
resistance), because the circuit reacts only to the real part of the input voltage, and the additional voltage drop on 
the radiation resistance Ur — jI 2 m Rrad is purely imaginary. 

An example topology of an electric dipole particle with the nonlinear active load is shown in Fig. 5, top, and the 
same for the magnetic dipole particle in Fig. 5, bottom (notice the difference in phases of the local oscillator signals 
in both schematics). We would like to stress that the operational amplifier seen in this schematic works with a DC 
signal. Thus, its role is not to amplify the evanescent fields of an incoming wave, but just to provide the necessary 
reactance of the loading circuit in order to tune the structure to a resonance. In turn, the evanescent modes in this 
structure are enhanced because of this resonance. 



V. CONCLUSIONS 

In this paper, the concept of perfect lensing with a pair of phase conjugating surfaces introduced by us earlier, 8 i.e., 
a possibility to achieve optical resolution well below the wavelength limit without using double-negative materials, 
has been further developed. Working as a planar lens, a pair of phase conjugating sheets is able to focus propagating 
modes of a source due to the negative refraction at the interfaces and, in the same time, enhance the evanescent modes 
due to surface plasmon-polariton resonances, i.e., it provides sub- wavelength resolution imaging indistinguishable from 
the perfect lens proposed by Pendry, while not suffering from its known drawbacks. 
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We have investigated in detail the physics of operation of nonlinear sheets with the boundary conditions of the form 
E t (w)|i = E t (w)*| 2 , H t (<x>)|i = H t (w)*|2 and have demonstrated that they are, in principle, physically realizable with 
devices imposing the necessary relations between the fields and the equivalent electric and magnetic surface currents 
at the phase conjugating boundary. Namely, we have shown that the mentioned surface currents must form Huygcns 
sources that radiate towards a given side of the boundary, negating the fields incident from the other side and creating 
complex-conjugated fields in the corresponding half space. 

As a possible realization of such sheets, we have proposed and considered in this paper arrays of nonlinear and 
nonreciprocal bi- anisotropic inclusions reacting differently to the propagating and evanescent plane waves. At mi- 
crowaves, the considered design makes use of balanced modulators (a type of mixers) to provide for the required 
nonlincarity and nonreciprocity of the circuit. At optics, a design utilizing similar principles may become feasible in 
future as the field of optical nanocircuits develops further. 

As a final note we would like to mention that in such arrays the enhancement of the evanescent waves is due to 
a high-quality surface mode resonance, as in the grids of passive resonant inclusions considered in Ref. 3. This is in 
contrast to Ref. 10 where the phase conjugating surface must itself parametrically amplify the fields which requires 
an unphysically high conversion efficiency. 
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